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in terms of the direct correlation function c(r, r'), the equilibrium density profile p (r) and an 
operator which relates density to surface configurations. Neither the concept of an effective interface 
Hamiltonian nor the difference in pressure is needed to determine the general structure of the height- 
height correlations or 7(q), respectively. This result generalizes the Mecke/Dietrich surface tension 
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Thermally excited capillary waves on the fluid inter- 
face between two phases have drawn the attention of 
many theoretical [ D, [ 3 , [3, B> H 0] and experimental 
BE HI Hi Hi USEE, EH physicists for more than 
hundred years. Two approaches have been developed in 
the last decades in order to resolve the structural prop- 
erties of the transition region and its thermally driven 
fluctuations: Density functional theory (DFT) is applied 
successfully in order to describe the dependence of inter- 
facial properties on, e.g., temperature or intermolecular 
forces, respectively (see, e.g., [13] )■ In particular, the 
spatial dependence of the one-particle equilibrium den- 
sity /3o(r), which exhibits a smooth transition between 
two coexisting phases, can be calculated within a DFT 
framework. 

On the other hand, capillary- wave theory (CWT), see 
Ref. analyzes the fluctuations of a infinitely thin and 
planar liquid-vapor interface parallel to the xy-plane 
M 2 . Local deviations from the planar surface config- 
uration are represented by a random field u(R), with 
R = (x,y) € R 2 , which is considered to be statistically 
independent of the particle distribution. Consequently, 
the height-height correlation function (u(R)u(R'))cwt 
depends parametrically on the surface tension 7 but does 
not include any particle interactions explicitly. 

In 1999, Mecke and Dietrich suggested a combination 
of DFT and CWT. They use an isodensity criterion, i.e., 
/o(R, z = u(R)J = const., in order to define the inter- 
face u(R) implicitly [ji]. By means of a grand canonical 
density functional Q [p] , they define an effective interface 
Hamiltonian 7i[u] := fl[p u ] — £l[po] as the difference be- 
tween two free energies, associated to two equilibrium 
surface states. 0,[p o ] is the free energy of a planar inter- 
face parallel to the xy-plane, while £l[p u ] is associated to 
a fixed but spatially varying surface it(R). An expansion 
of Q[p u ] near Q[p Q ] (or p u (r) at p {r), respectively) al- 
lows for an explicit expression of TL[u] in terms of w(R). 
By taking curvature corrections of the fluctuating inter- 
face into account, the resulting height correlation func- 



tion is governed by a wavelength dependent surface ten- 
sion 7md(<z)- While the included van der Waals attrac- 
tions induce a decrease of 7md(<?), its increase at large 
values of q is considered as an indication for a bending 
rigidity for liquid surfaces. 7md(<?) has been confirmed in 
x-ray studies for several liquids with different accuracy 
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Although the combination of DFT and CWT has been 
generalized to binary mixtures Q , the results depend on 
the particular choice of the density functional Cl[p] and 
on the validity of the curvature corrections. There exist 
also difficulties in analyzing the exp erimental scattering 
data in order to obtain 7md(<?) uM- In addition, recent 
numerical studies even cast doubt on 7md 
Ref. [3] the authors argue that 7md(<z) must decrease for 
large values of q if derived from a DFT as in Ref. [f|. 
Thus, alternative theoretical methods are mandatory in 
order to elucidate the relation between density correla- 
tions and interface correlations. 

Here the height-height correlation function for a fluctu- 
ating interface between two coexisting phases is derived 
and expressed by means of general equilibrium proper- 
ties of the corresponding density-density correlations. In 
principle, this approach is neither restricted to planar 
systems nor to liquid- vapor interfaces. An explicit ex- 
pression for a wavelength-dependent surface tension 7(q) 
can be identified, which depends on the direct correlation 
function c(r,r') (which is defined via the non-ideal gas 
contributions of the inverse density covariance function), 
the equilibrium density profile p (r) and, an operator re- 
lating density and surface configurations. The expression 
for 7(q) generalizes in particular 7m d (9) and its deriva- 
tion does not rely on a free energy expansion, which is 
one of the main criticism of Ref. [19( . 

We consider an equilibrium state within the grand 
canonical ensemble. Ensemble averages are denoted 
as („.}o. Using the local microscopic density p(r) :— 
YliLi <K r ~ r we write jOo(r) := (p(r))o for the spa- 
tially dependent one-particle equilibrium density. We 
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assume, that the system exhibits a single interface be- 
tween two coexisting phases. Its mathematical speci- 
fication is given below. Density fluctuations are char- 
acterized by the density covariance function G(r, r') := 
( (p(r) — p (r)) (p(r') — p (r')) ) . Since fluctuations of 
the interface are not independent of density fluctuations, 
the height-height correlation function must be related to 
G(r,r'). In the following, we study this relationship. 

It follows from probability theory, that a random 
field X(r) G M exists, which exhibits the same mean 
value and correlations as p(r), i.e., X(r) = p (r) and 

[X(r)-p (r)] [X(r') -Po(r')] = G(r, r'), where the av- 
erage "777 is taken with respect to a suitable probability 
measure for X 20] . In general, the system under consid- 
eration is inhomogeneous due to the presence of an in- 
terface and thus G(r, r') does not depend only on r - r'. 
Since the first and the second moment of X(r) are given 
by p (r) and the G(r,r'), X(r) can be considered as, 
e.g., an inhomogeneous Gaussian process. By construc- 
tion, the local microscopic density p(r) and the random 
field X(r) can be distinguished only by measurements or 
calculations of their higher correlation functions. In par- 
ticular, the structure factor is not sufficient to determine 
the difference between X(r) and p(r). Thus, X(r) is not 
an equilibrium density but mimics the first and the sec- 
ond equilibrium moment of p(r). We call X(r) a random 
density. If, in addition, X(r) is generated by another 
random field which mimics the interface configurations, 
the correlations of this underlying process are related to 
G(r, r'). In the following, this procedure and the relation 
between the different correlation functions is discussed. 

Although we will focus later on planar interfaces, we 
use first a slightly more abstract notation in order to 
isolate the core definitions and ideas from additional as- 
sumption about the system like its symmetry. 

Corresponding to the equilibrium density /0 o (r) we in- 
troduce an equilibrium phase boundary s (R) G R 3 , 
which depends on two parameters R G M C R 2 for 
a suitable subset M of R 2 . The interface s (R) is 
considered to be an isodensity contour of p (r), i.e., 
/ o (s o (R)J = p c = const, for all R G M. Introducing the 
random field u(R) G R, a fluctuating interface s(R) G R 3 
is considered to be a random field which results from 
s (R) by local normal displacements (see Fig. [I}, 

s(R) = s (R) +u(R)n G (R) , (1) 

where n (R) is the unit normal vector of the surface 
s (R). For instance, we expect n (x, y) = e z for a liquid- 
vapor interface in a homogeneous field acting along the 
z-axis e z , while for a spherical droplet it is n D ($, ip) = 
e r (i?, Lp), where e r the radial unit vector. 

A relation between the random density X(r) and the 
random surface s(R) establishes a relation between the 
density correlations G(r, r') and the correlations of w(R). 




Figure 1: The random surface s is generated by random dis- 
placements u(R) of s (R) along its normal vector n (R) 
(see Eq. |T}). For simplicity, the surface patch s is drawn 
spherically. For each spatial point r we write d = d(r) 
for its (normal) distance to the equilibrium interface s (see 
Eq. ((3])). Random deviations of the equilibrium density p (r) 
are caused by the random displacements u(R) of the interface 
as well as by a local change in volume (see Eqs. Q and ([6])). 

For this purpose, we require the following conditions: 

A(s(R)) = po(s (R)) = Pc (2a) 
s = s X(r) = p Q (r) for all r , (2b) 

where s = s means the identity for all R G M. The 
condition Eq. (f2"a|) generalizes the concept of an isoden- 
sity contour to X(r). The condition (|2b[) states that X(r) 
assumes its mean value if s is identical to s . Thus, we 
consider those X(r) which result from p (r) due to a 
change in the isodensity surface. 

The above-mentioned conditions are not sufficient to 
determine a unique relation between X(r) and s(R). 
But here we show in which way different realizations 
of Eqs. lead to different height correlations for it(R) 
exploiting only the above-mentioned condition that the 
auto-correlation of X(r) — p (r) is given by G(r, r'). 

We assume, that each point r can be represented un- 
ambiguously within the normal coordinate system of the 
equilibrium surface s , i.e., there exist d(r) = d G K and 
R(r) = R G M so that (see Fig.HJ (Hf 

r = s o (R) + dn (R) . (3) 

A simple choice for X(r) which fulfills Eqs. ((2|) is 

X(t) = p o (r + s (R)-s(R)) , (4) 

where R = R(r) (see Eq. ©) By Eq. Q, X(r) is 

generated by random displacements of the equilibrium 
density p (r) along the normal vector n in whose direc- 
tion r lies. It is worth to note, that Eq. ((4]) is not a result 
but rather a choice for the parametrization of X(r) based 
on Eqs.©. Although Eq. ((4|) appears as the standard 
approach in order to describe capillary wave fluctuations 
(of planar interfaces), the specification of the conditions 
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Eqs. ([2]) is an attempt to justify this approach on a more 
profound principle. In particular, it does not rely on the 
thermodynamic conditions as, e.g., the Gibbs dividing 
surface concept. Here, we are interested the generaliza- 
tion of Eq. ^ which are consistent with that principle. 

In general, a change in density at point r is obtained 
due to a change of the infinitesimally small volume which 
contains the point r. This change of the local volume at 
point r might occur from a change of the equilibrium 
interface even if the corresponding point on s D is not 
shifted along n . For instance, a change from a spherical 
droplet to an ellipsoid changes the local volume at a point 
r without moving the corresponding point on s . Thus, 
a change in curvature is one of these mechanisms. Obvi- 
ously, such changes in local volume depends on the dis- 
tance d between r and s (see Eq. (j3|) and Fig.[lJ. Due to 
the isodensity condition this effect vanishes for d = u(R), 
i.e., if r lies on s(R). Equally, if the amplitude of the 
interface fluctuations are sufficiently small, the local vol- 
ume within the bulk phases does not change. Thus, this 
effect vanishes for d » £, where £ is the effective width 
of the interface, i.e., the width of the transition region in 
which p (r) differs considerably from the bulk densities. 
For step- like profiles it is £ = 0. In such cases we do not 
expect such a change in local volume due to other mech- 
anisms. A change in local volume that affects the local 
density p (r) can be associated with a local compressibil- 
ity. Therefore, the following generalizations of Eq. ^ are 
intended to take into account the local compressibility. 

The incompleteness of the standard capillary wave 
ansatz Eq. ([4]) can be understood equally from the phys- 
ical picture, that density fluctuations are present in each 
spatial direction and not only along the normal direction 
of the interface. This lack in description could be ad- 
justed by taking tangential fluctuations (parallel to the 
interface) into account. On the other hand, tangential 
density fluctuations affect the normal density fluctua- 
tions due to the isodensity condition of the interface. 
Therefore we expect an additional term in Eq. (fj[|) which 
projects tangential density fluctuation on normal den- 
sity fluctuations, similar as it is done in the projector 
operator formalism. This projector includes the ratio of 
correlation functions associated with tangential and nor- 
mal fluctuations, respectively. Furthermore, it depends 
in general on the distance to the interface and on the 
wavelength of the fluctuation: For large distances the 
tangential density fluctuations are not expected to affect 
the normal density fluctuation which implies that the 
projected contribution disappears independently of the 
wavelength of the fluctuation. Close to the interface, a 
long- wavelength tangential density fluctuation will barely 
have an influence on the local interfacial structure due to 
the isodensity condition, i.e., such tangential fluctuations 
shall be rather aligned by the global structure of the in- 
terface. Thus, the projector is expected to vanish in that 
case. For the same reason, the local interfacial struc- 



ture depends on short-wavelength tangential fluctuation 
close to the interface. For instance, the interface might be 
bend locally due to a tangential density fluctuation above 
or below the interface. At this point the interpretation 
in terms of tangential density fluctuations becomes simi- 
lar to the above-mentioned local compressibility picture. 
Therefore, the interfacial width £ can be considered also 
as the range of the coupling between these two types of 
fluctuations. 

After these pictorial remarks we introduce the in- 
dependent linear operator L(r) which acts on the surface 
s(R) (and s (R)) via 

L(r)s(R) := f L(r; R, R') s(R') dR' , (5) 

J M 

with a 3 x 3-matrix L(r; R, R') as integral kernel. Since 
L(r) is intended to model the effect of local compressibil- 
ity we expect L(r)s (R) = for all R 6 M and all r 
due to the isodensity condition. Thus, it is L(r)s(R) = 
L(r) (u(R)n G (R)). Equally we assume L(0) = 0. For 
this reason, a more general choice for X(r) which fulfills 
the conditions in Eqs. (J2J) is 

X(r) = Po (r - u(R)n (R) + L(r - s(R))s(R)) , (6) 

where R = R(r) (see Eq. ([3])). As mentioned above, 
this study elucidates the general influence of L(r) on 
the height correlations. We will not derive its partic- 
ular form. From the physical meaning discussed above 
in terms of tangential density fluctuations, it becomes 
clear, that the integral kernel L(r; R, R') should reflect 
the presence of an interfacial symmetry of the equilibrium 
interface s . We will demonstrate this principle below in 
the case of planar interfaces. 

An expansion of X(r) given by Eq. ^ with respect to 
it(R) up to linear order gives 

X(r) - Po (r) w - Vpo(r) [l 3 - L(dn (R))] u(R)n (R) , 

(7) 

where I3 means the 3x3 unit matrix, d = d(r) is the 
distance of r to the surface s and R = R(r) (see Eq. Q 
and Fig.[T|). L(dn (R))u(R)n (R) can be decomposed 
into normal and tangential contributions with respect to 
s (R) at each R which allows for a general derivation of 
the height correlations. Here, we demonstrate the main 
idea for planar interfaces. More general cases like the 
spherical or the cylindrical interface differ by their coor- 
dinate system which makes primarily the notation more 
complex. In addition, closed interfaces (e.g. droplets) are 
subject to an additional volume constraint for stability 
reasons; these exclude certain types of fluctuations. The 
crucial point in all cases is the symmetry of the equilib- 
rium interface which is reflected by the correlation func- 
tion G(r, r') and the integral kernel L(r; R, R'). 

In the following discussion of height correlations of 
infinite planar interfaces (R = (i?i,i?2) £ M = R 2 ), 
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the equilibrium surface s (R) = R\ e x + R% e y + z Q e z , 
with a constant z a , lies parallel to (x, y)-plane so that 
n Q (R) = e 2 . The equilibrium density p (r) = p (z) de- 
pends only on the normal distance z = z D + d(r) to the 
surface s (see Eq. $5$ and Fig.Q}. Consequently, the 
density correlation function G(r,r') = G(z,z',R — R') 
is homogeneous with respect to the lateral coordinates. 
Equally, the integral kernel possesses the same homo- 
geneity, i.e., L(dn ;R, R') = L(z — c;R — R'). Since 
p (z) depends only on z, it is d x p — d y p = and 
Eq. ([7]) becomes 

X(z, R) - p (z) « - d zPo (z) [1 - l(z - z )} u(R) , (8) 

where L(z - z D )u(R) = / R2 L(z - z ;TL - R) u(R') dR 
with L(0) = is a linear operator which depends on z — z 
and which acts on the random variable w(R) (see also 
Eq. ([5])). The correlation function for X(z,R) — p (z) 
is given by G(z,z',R — R'). By taking the Fourier- 
transformation with respect to the lateral coordinates, 
we obtain from Eq. (JSj> 



27r5(q + q')G(z,z',q) 
d z p (z)d z p (z') x 



(9) 



[1 - L(z - z Q ; q)][l - L(z' - z ; q')] fi(q) fi(q') . 

The interpretation of L as a projector of tangential 
density fluctuations on normal density fluctuations im- 
plies, that L includes the ratio of Fourier components 
of the transverse auto-correlations and the height auto- 
correlation ?i(q)u(q') . Therefore and from the discus- 
sion above it becomes clear, that, in particular for larger 
q values, G is reasonably governed by the transverse cor- 
relations and not only by the height correlations as in the 
case L = 0. 

G(z,z',q) is a positive definite function and symmet- 
ric in z and z'. This allows for an expansion of Eq. Q 
in terms of eigenfunctions of G(z,z',q), which leads to 
a generalized version of Wertheim's eigenfunction anal- 
ysis of the correlations in a planar liquid-gas interface 
for small values of |q| 24]. The need of such a gen- 
eralized approach has been already concluded from nu- 
merical investigations of the eigenfunction ansatz [25| . 
Here, we first ask in which way the term L(z — z ;q) 
affects the general form of u(q)u(q') which will be an- 
swered below without using an expansion in eigenfunc- 
tions. For low q values, the result becomes equal to the 
expression from Wertheim's approach. It is important 
to bear in mind, that G(z, z', q) results from a thermo- 
dynamic average procedure while the height correlation 
function u(q)w(q') stems from an unknown probability 
measure. Therefore, Eq. §§§ is in particular a manifes- 
tation of the required thermodynamic consistency condi- 
tion mentioned above. 

The inverse density correlation function G _1 (r, r') is 
defined by J dr' G(r, ^G' 1 ^, r") = 6{r - r"). This 



implies from Eq. ([9]) for the height correlation function 
u(q)u(q') = 2 7 r ( 5(q + q')G(q) (10a) 
G(q) := 

[1 - L(z - z ; - 



dzdz' d z p {z)d z > p {z') x (10b) 
q)}[l-L(z' -z ;q)}G-\z,z',q) 



G _1 is typically written as G _1 (r, r') = 6(r— r')/p (r) — 
c(r , r'), where c(r, r') is called direct correlation function 
[261 ]. In Fourier space, we have for the planar interfaces 
G _1 (z, z', q) = S(z — z')/p (z) — c(z, z', q). Furthermore, 
the equilibrium density p (z) fulfills the generalized baro- 
metric law, i.e., the equation \n\^ h p (z) — c^(z) + 
/3V ext (z) = [Sfj,, where A t h is the thermal de-Broglie wave- 
length, c^ x \z) means the effective one-particle potential, 
V ext (z) is an external potential and, p means the chem- 
ical potential Q. By taking the derivative of the equi- 
librium condition for p Q (z) and bearing in mind the re- 
lation VcW(r) = J dr'c(r,r')VVo(r') (see Ref. 0) we 
can rewrite the ideal gas contribution in G _1 . This gives 



G- 1 (z,z',q) = -/3 



6(z-z')d z V ext (z) 



d z p (z) 



c(z,z',q) (11) 



5(z - z') 



dz" c{z,z",0)d z p o (z") 



d z p {z) 

By combining Eqs. (flU)) and (fTTj) we obtain 

27T^(q + q') 



P u(q)tt(q') 



u ext (q)+7?(q) 



where 

« ext (q) 

and 

?y(q) := 



dzd z V ext (z)d z p (z)\l-L(z 



dzdz' d z p (z)d z > p (z') x 



(12a) 



z ;q)\ 
(12b) 



(12c) 



(c(z,z',0) - c(z,z',q)) |l - L(z- z ;q)\ 



+ - c(z, z', q) \L(z - z 



q) - L(z' - 2o;q)| 



By construction, the height-height correlations given by 
Eqs. (H2) are thermodynamically consistent with the den- 
sity covariance function G(r, r'). From the derivation 
above the q-dependence of r)(q) comes in via the direct 
correlation function c(r, r') and the kernel L(z; q) of the 
linear operator L(z). While c(r,r') represents the parti- 
cle interactions, L(z) takes into account the change of the 
equilibrium density p (z) due to a local volume change 
at point z, i.e., the local compressibility, or the ratio of 
tangential density correlations and normal density corre- 
lations, respectively. 

The derivation of the explicit expression for 
u(q)u(q'), i.e., from Eq. © to Eqs. (fl"2"|) , is based 
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on the definition for G _1 (z,z',q) and the generalized 
barometric law. Therefore, any choice for L leads to a 
height correlation function Eqs. p"2"|) which is consistent 
with the density covariance G(r, r'). Putting Eq. Q12ap 
back into Eq. © it follows 



(3G(z,z',q) = dp (z)dp (z') x 



[l-L(z-z ;q)}[l-L(z 



(13) 

-q)] 



« ext (q) + ^(q) 

where the functional dependence of w ext (q) and rj(q) on 
L(z ~ z ; q) is given by Eq. (|12b[) and p2cp . respectively. 
Consequently, for a given G(z, z', q) the Eq. p3|) is a 
definition for L(z — z ;q) (possibly limited by the lin- 
earization made in Eq. (j8])). For instance, using the total 
correlation function h(z, z', R — R'), it is G(z,z',q) = 

(5(z — z') j o (z) + j0 (^)jOo(^')' 1 ( 2; j z '7 q)- From the Ornstein- 
Zcrnikc relation for inhomogeneous systems combined 
with some closure relation, an approximate scheme for 
h(z,z',q) can be applied in order to get a functional 
equation for L{z — z ; q). In reverse, a given L implies a 
particular form of G(z,z',q), where the physical mean- 
ing of L (as discussed above Eq. might serve as a 
guidance. The corresponding height correlation function 
is in any case consistently given by the Eqs. (fT2"]) . That 
means in particular for larger q values, that a sole investi- 
gation of u(q) w(q') , i.e., a sole examination of Eqs. (fT2| . 
is not sufficient in order to check the reliability of a given 
L(z-z ;q) [27|. 

The only approximations in the realization of Eqs. ([2]) 
are those in Eq. © or ((7J) , respectively. Surprisingly, the 
concept of an effective interface Hamiltonian or an ex- 
pansion of free energies, respectively, is not needed in 
order to derive the general structure of height correla- 
tions which are consistent with the underlying density 
correlations. Nevertheless, the Eqs. (fl2|) and in particu- 
lar Eq. (|12c| can be compared to former expression for 
the height correlations based on such principles. To do 
so, it is convenient to consider j(q) :— r)(\q\)/q 2 which is 
referred to in literature as a wavelength dependent sur- 
face energy density for isotropic interfaces. 

The macroscopic surface tension 7(0) has been derived 
in Ref. [28] and follows from Eq. (TT2cD, if L{z- q) -> for 
<7 ^ 0. For step-like profiles z p o (z) ~ 5(z — z ) the 
the L-terms drop out due to L(0, q) = 0. In that case, 
Eq. (|12c|) reduces to the result derived in Ref. [4| if the 
same expression for direct correlation function provided 
by the density functional in [1] is used (see Eq. (4.6) in 
Ref. H). Both cases reflect some properties of L which 
we concluded from the general discussion above, i.e., the 
long-wavelength limit and the £ = case. 

The Eqs. (|12p include also the energy density for sur- 
face excitations derived by Mecke and Dietrich in Gaus- 
sian approximation 1. To see that, we first note, 
that the particular density functional used in Q im- 
plies an explicit expression G^q(z, z', q) for the inverse 



density covariance function and thus for the direct cor- 
relation function CMo(z,z',q) = S ^ z - ) - G^ D (z,z' ,q) 
[29I ] . In Ref. the external potential is the homoge- 
neous gravity potential V^q(z) = mgz and the equi- 
librium interface is located at z = (see Eq. (2.4) 
in 0). By repeating the derivation of the Eqs. (TT2")) 
from Eq. |TUD, one applies Eq. $11]) with V ext (z) = 
V^(z) and c(z,z',q) = cmd(^, z' ,q) except for the 
term d z p (z)d Z ' p (z') L*(z:q)L(z';q)G~ 1 (z,z / ,q) which 
appears in Eq. (1 10[) . Therein, one uses G[^q(z, z', ^) in- 
stead of Eq. pip (which leads to the n — k (q) contri- 
bution in Eq. (3.11) in [6[). By setting d z p (z) L(z; q) = 
— Q 2 Ph(z) (where pn is given by the Eqs. (3.27) and 
(3.31) in Ref. Q without further derivation) the result- 
ing expression for w(q) u(q') and n{q)/q 2 are equal to 
Eqs. (4.1), (4.2) and (3.11), respectively, in Ref. Q (see 
also Eqs. (2.28) and (3.10) in 6]). From the explicit ex- 
pression for ph(z) and p (z) assumed in Q we obtain 
iM D (z,q) = sinh(^)q 2 [30]. Eq.'CE]) implies, 

that Lmd(z, q) combined with cmd(z, z' , q) provide rather 
a particular model for the density covariance function 
G(z, z' , q) than a model for the height correlations, only. 
Therefore, in order to test the reliability of the resulting 
expressions in computer simulations it might be easier to 
check Eq. (|T3j) instead of Eqs. (fl"2"]) via arguable numer- 
ical procedures to define the position of the fluctuating 
interface. 

The form of Lfj\o(z,q) implies, that the operator 
Lmd( z ) is not bounded for z / 0, since Lmd(z) ~ V 2 so 
that the kernel Lmd(z, R— R') involves the ^-distribution 
and its second derivative in x- and ^/-direction. From 
Eq. ([6]) follows that X(r) E [ptiPo] where pf := p (z — > 
±00) denote the bulk densities. Therefore, a realiza- 
tion w(R) of the interface with u(R p ) = |Vu(R p )| = 
and |V 2 u(R p )| = 00 for a particular point R p induces 
X(z ^ 0,R P ) = pf [Hj]. On the other hand, the lin- 
earization of Eq. (J6j> (Eq. ([8]) for planar interfaces) with 
L = I_md allows for X(r) G [— 00, 00] even if it(R) is 
bounded, which seems not reasonable from the physical 
point of view. Therefore, as long as one would like to lin- 
earize Eq. ©, one has to limit the applicability of Lmd(-z) 
to a particular set of interface configurations whose sec- 
ond derivatives are also bounded or equivalcntly, one has 
to limit the g-range of the kernel LMD(z,q) [32] - 

The limitation of the set of interface configurations 
{w(R)} restricts the set of modeled density configura- 
tions. If {X(r)} denotes the set of all density configu- 
rations, an optimal sampling of {X(r)} would capture 
the set of all relevant density configurations {X(r)} re i C 
{X(r)} which, in our case, contribute to the density co- 
variance function G(r,r'). The representation of X(r) 
by the interfacial field w(R) leads to a set of density con- 
figurations {X(r)}\f which gives in general not the opti- 
mal sampling of {X(r)} such that {X(r)}\f C {X(r)} re \. 
Therefore one should allow for the largest set of inter- 



6 



face configurations {it(R)} (whose Fourier transform is 
continuous, for instance) in order to exhaust {X(r)} re | 
in an optimal way, which corresponds to an approxima- 
tion of G(r, r') by the height correlation function via 
Eq. (|13[) in an optimal way. As an example, one may 
allow for all interface configurations u(R) with ||u||i := 
L 2 |u(R)| dR < oo and in order to use Eq. J5]), one may 
require \L(z, R)| < oo for all R £ M 2 so that L(z) be- 
comes a bounded operator for every \z\ < oo [33j. 

The operator Lmd(z 7^ 0) ~ V 2 is associated to the lo- 
cal curvature of the interface u(R), since it results from a 
truncated curvature expansion of the corresponding den- 
sity configuration (see Eq. (2.17) in [6]). Consequently, 
the coefficient of the resulting g 2 -increase in 7m d (9) has 
been interpreted as the bending rigidity of the liquid in- 
terface. The above considerations about the boundedness 
of L(z) then suggest, that the influnce of the local cur- 
vature of the interface and the effect of bending rigidity 
can be probably not continued to arbitrary small scales 
although the random density X(r) as well as the random 
interface w(R) are defined on all length scales. Neverthe- 
less, since L(z)u(R) is written as a convolution integral 
(see Eq. ((8]) below) and L(z) is bounded, we thus may 
imagine that L(z) picks up the non-local effects of the 
interface u(R) on the density configuration X(r). In a 
similar manner, the influence of non-locality has been 
mentioned also in the framework of short-ranged wetting 

In our treatment of interface fluctuations, the dis- 
torted interface s(R) is considered to be generated by 
local random displacements w(R) along the normal vec- 
tor n (R) of the equilibrium interface s (R) (see Eq. JT]) 
and Fig. (p}). By linking the random displacements to 
a change in density via Eq. ([6]), the density correlations 
G(r, r') govern the interface correlations u(R)u(R') . 
The situation is in a sense reverse to those in stochas- 
tic differential equations, where the mean value and the 
correlations of added noise terms (which correspond to 
m(R) here) are specified in order to study its influence on 
a variable of interest (which is p(r) in our case). In gen- 
eral, the noise contribution within a single realization of 
the process can not be determined unambiguously. Simi- 
larly, the difficulty to obtain u(R)it(R') from numerical 
simulations arises from the problem to identify the ran- 
dom interface, i.e., the realization of the related (lower 
dimensional) stochastic process, in a particle configura- 



tion [22] , |23| . Consequently, any recipe that determines 
the random interface in computer simulations works ac- 
curately with a certain probability but can not be exact. 
As pointed out in Ref. [l9(, the precision of such proce- 
dures is enhanced the more microscopic information of 
the system is included, i.e., the more n-particle correla- 
tions are taken into account to define the random inter- 
face numerically. 

At this point we come back to the criticism of 
Tarazona, Checa, and Chacon [lj| concerning the 



Mecke/Dietrich approach. In order to obtain the same 
expression for the height correlation function as pub- 
lished in @ we have neither used an effective interface 
Hamiltonian nor a restricted variational principle as as- 
sumed in Q. Formally, the random density X(r) cor- 
responds to Pf(r) in Ref. @ (see Eq. (2.6) in [||) but 
the crucial difference consists in the assumption, that 
Pf(r) in [f| minimizes the given density functional with 
an additional isodensity condition (see Eq. (2.5) in 0]). 
While the isodensity condition for pj in 6] is similar 
to the isodensity condition Eq. (f2"a[) for X(r), no addi- 
tional equation resulting from a minimization procedure 
is required for X(r), i.e., the realization of Eqs. ((2]) are 
those in Eq. (J6J) or ([7j , respectively. In order to extract 
the height correlation function or 7md(?), respectively, 
from computer simulations, the authors of Ref. [l9| have 
taken this minimization condition for pf in Q seriously 
into account. Their results differ considerably from the 
predicted 7m d (9) due to general features of density func- 
tionals as explained in Ref. flit ]. Since pf fulfills a min- 
imization condition for a density functional, it includes 
capillary waves on small wavelengths which can not be 
separated as height fluctuations (see Fig. 2 in [l9j]). While 
the numerical analysis shows, that the minimization con- 
dition leads to a different 7md(?) than predicted, our ap- 
proach shows, that 7md(<?) can be derived indeed without 
an additional (minimization) condition. In other words, 
the results in [19( do not necessarily imply that 7md(<?) 
is structurally incorrect because no density functional 
and related minimization procedures are needed in or- 
der to derive 7md(?)- The only dubious quantity that 
remains is L(z — z ,q) or Lmo{z, q), respectively, and 
from the derivation above it becomes clear, that a model 
for L(z — z a ,q) is a different task, that can not be solved 
within a capillary wave theory. In particular, our deriva- 
tion implies that L{z — z Ql q) or Lmo(z, q) should be dis- 
cussed rather in terms of Eq. ()13|) than only in terms of 
Eqs. (HU). This suggestion if fully consistent with the con- 
clusion in [lljj . that the capillary wave problem can not 
be solved (numerically) on short wavelengths by taking 
into account only the one particle distribution. 

Above we have discussed in detail the g-dependence 
of L(z,q) or Lmd(-z,<7), respectively. Finally we briefly 
discuss a recently published study on planar colloid- 
polymer interfaces by Blokhuis, Kuipers and Vink [35I ] 
as an application of Eq. (|13|) (with z = 0) in order 
to show that also the z-dependence of L(z, q) plays 
a crucial role. In Ref. [35j a particular model for 
the so-called surface density- density correlation function 
S(q) := (Ap)- 2 G(z, z', q)dzdz' (see Eq. (1) in 

Ref. [35| where Ap := p~ — p + is the difference of the 
coexisting bulk densities) is proposed in order to ex- 
plain the numerical data. SBKv(q) m contains an 
interfacial contribution Sp.h(o) and a bulk contribution 
S b (q) (see Eq. (11) in [35j). Since S'b(q) is modeled in- 
dependently from Shh{q) (the bulk density fluctuations 
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are considered to be uncorrelated from the interface fluc- 
tuations, see Eq. (10) in (35[) and only Shh(o) is used 
to adjust the data fit (Eqs. (14) and (15) in [35|), we 
focus here only on SV^q). In order to compare S(q) 
resulting from Eq. (fT3")) to S^/^q), we adpot v ext = 
and use also the Helfrich form rj(q)/q 2 = a + q 2 n (see 
Eqs. (8) and (9) in [35|), where a means the macroscopic 
surface tension and n > means the (phenomenologi- 
cal) bending rigidity. From Eq. (4) in [35[ we identify 
d z p (z)L B Kv(z, q) = \pi{z)q 2 , where pi(z) is given by 
Eq. (26) in j3f|. Interestingly, the kernels I/MD(z,q) and 
LBVK(z,q) show the same g-dependence, but the weight 
functions Ph(z) and pi(z), respectively, are different. In 
particular, it is p\{z — 0) ^ so that the isodensity 
condition Eq. (|2a|) is not fulfilled. By evaluating the in- 
tegral / d z p (z) (l — L(z, q)) dz, we then obtain from 
Eq. (TT3D 



S(q) 



(1 + Cq 



2\2 



+ f q 2 ) 



(14) 



with C = C H f 2 (MD) or C = (BKV), respectively 
(seejpl and Eq. (6) in 35]). The best fit to the data 
in J35| for in the entire q-range results from Eq. (15) in 
[Ijwhich implies S hh {q) = iPaq 2 )" 1 (1 + f q 2 ) with 
k = — kbkv > because kbkv is found to be negative 
(Table I in 35]). This form for 5^h(q) is obtained from 
Eq. (fT4|) for K = CaiiC>0 which is indeed the case 
for the Lmd(z,<i) but not for Lbvk(z, q). Since Eq. JT 
expresses the consistency between the height correlation 
function and the density covariance function, this re- 
sult indicates that the kernel Lmd(z, q) seems to be even 
more consistent than Lbvk(z, q) with both, the assumed 
Helfrich form for r](q)/q 2 and the numerical results for 
Shhiq) for the hole q-range in [35]. In particular, the 
negative bending rigitity kbkv in [3a| appears rather ar- 
tificially from the truncated Taylor expansion of Shh(q) 
in Eq. (14) in [2E], which is assumed to be applicable for 
all q. Another main result of (Hf is, that kbkv vanishes 
for T — > T c . Using k = CrO'^ 2 from the relation above, 
it then follows that Ch vanishes with T — > T c . This is 
in qualitative agreement with Cjj(T) ~ £ _S (T) for s > 
which has been concluded from general considerations 
about the expected quantitative influence of Lmo(z, q) at 
higher temperatures (see comments on Eq. (55) in [7J] and 
discussion on Eq. (4.72) in 13610 . but without determining 
the exponent s rigorously [37j. Thus, although our ap- 
proach does not provide a particular theory for L(z,q), 
it might serve as a tool in order to verify the consistency 
of various models with the underlying density covariance 
function. 



X-ray experiments indicate L(z;q) ~ q 2 for isotropic 



interfaces 0, ES EH El 0] • But a separation of the 
height correlations from the scattering data corresponds 
to the Eq. (TT3")) . As discussed above, any choice for L 
leads to a height correlation function which is consistent 



with G. This makes probably clear uncertainty of such 
procedures at larger q values. Although a generalization 
of Eqs. (fT2| to multi-component liquid mixtures as well 
as for spherical interfaces can be performed, the more 
challenging task is to identify the correlation ratios or 
the local compressibility, respectively, hidden in L. 
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and P. Tarazona. This work has been supported 
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